Decoherence and maximal violations of the Leggett-Garg inequality 
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We consider maximal violations of the Leggett-Garg inequality, obtained by maximising over all 
possible measurement operators, in relation to non-unitary aspects of the system dynamics. We 
model the action of an environment on a qubit in terms of generic quantum channels and relate 
the maximal value of the Leggett-Garg correlator to the channel parameters. We focus on unital 
channels, and hence on decoherence. In certain important cases, exact relations between the channel 
parameters and maximal violations can be found. Moreover, we demonstrate the existence of distinct 
thresholds for the channel parameters, below which no violation of the Leggett-Garg inequality can 



PACS numbers: 03.65.Ud, 03.65.Ta, 

I. INTRODUCTION 

The Leggett-Garg inequality (LGI) pj has been the 
focus of much recent attention, both in theory [H-El and 
in experiment fL2Hl9j . The inequality (or its simplest 
variant) reads 



K = C21 + C 32 - C31 < 1, 



(1) 



where C a p = (Q a (ta)Qp(t/3)) is the correlation function 
of the dichotomous variables Q a ,/3 = ±1 at times t = t a 
and t = tp. Under the condition that the measurements 
are performed non-invasively, violation of this inequal- 
ity implies the absence of a macroscopic real description 
of the system [l| . Due to the existence of superposition 
states, quantum mechanics is not a macroscopic real the- 
ory, and can thus violate Eq. (JXJ) . A two- level system 
undergoing Rabi oscillations, for example, can achieve a 
value of K = | for appropriate parameters. Violation of 
the LGI by a quantum system is related to, and is typi- 
cally taken as an indication of, quantum coherence in the 
system. 

It is the aim of this paper to make a quantitative link 
between coherence — or rather, decoherence, the process 
by which coherence is lost — and violations of LGI. The 
point-of-view we take is that, if we had a perfect quantum 
system and could control its evolution as we please, then 
we should be able to violate Eq. ([T]) maximally. On the 
other hand, a system which supports no coherence (over 
the relevant time-scales) should certainly not produce a 
violation. In between these two extremes lies a spectrum 
of partially coherent dynamics and the questions is what 
happens for such intermediate cases? What quantitative 
statements can we make concerning the degree and kind 
of decoherence and violations of LGI? 

Our approach to answering these questions involves of 
two elements: First, we describe the time evolution of 
the system between measurements in general terms by 
making use of the formalism of quantum channels [20j, 
l2lj |. Second, we shall study the maximal violations of 
Eq. ([I} obtained by maximising over all possible choices 
of the measurement operators Q a . This approach takes a 
significant part of the unitary component of the quantum 



evolution as given and allows us to connect if maX 7 the 
maximised value of the correlator from Eq. to the 
non-unitary aspects of the dynamics. 

The template for this approach is that taken with 
the (spatial) Bell's inequalities [H In the CHSH 

inequality [2^, for example, maximising over all possible 
detector angles establishes a direct relationship between 
the Bell correlator and the entanglement of the state un- 
der test (at least for pure state of two qubits) [25j . For 
the CHSH inequality, maximization over measurement 
angles connects the value of the Bell correlator with a 
property of the input state (entanglement). In contrast, 
maximization over measurement angles for the LGI re- 
veals the connection between the Leggett-Garg correla- 
tor, K, and a property of the dynamics (decoherence, 
dissipation) . 

Whilst our approach is general, we will focus here on 
a single qubit and environments whose effect on the sys- 
tem are Markovian. We shall also focus on situations in 
which the system undergoes decoherence but not relax- 
ation, i. e. the channels we will mostly consider arc unital 
channels J2(j. Non- unital channels are briefly discussed 
in Sec. ED 



II. LGI IN TERMS OF QUANTUM CHANNELS 

We begin by describing our generic model of the mea- 
surement of Eq. ([!]) for a qubit. 



A. Measurements 

We assume that the measurements used to build C a p 
are projections performed in an arbitrary basis. Specifi- 
cally, we parameterise the measurement of Q a as first a 
rotation of the system with the unitary 
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(2) 



followed by a projection on to a z cigenstates, and then 
by the back-rotation, q\. The measurement procedure 
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can the be described by the map p — > Q a [p\ acting on 
density matrix p. 



B. Quantum channels 

We will describe the time evolution between measure- 
ments in terms of the formalism of quantum channels 
[20I |2lj . The most general map of one density matrix 
onto another, p — > $ [p], can be written 

$ [p] = v*S [upv)\ v, (3) 

where u and v are unitary rotations and S is a supcrop- 
erator responsible for non-unitary evolution whose action 
can be specified as follows. Every density matrix can be 
written p — i (1 + w.cr) where u is the vector of Pauli 
matrices and w is a real, length-3 vector with |w| < 1. 
Representing p by length-4 vector \p)) = (1, w) the action 
of S can be expressed in terms of a matrix acting on this 
vector: S[p] <-» S\p)). The most general form of S is 
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(4) 



with bi and c, real. That the map must be stochastic (i.e. 
S[p] is a valid density matrix) implies constraints on the 
coefficients bi and Cj. 

The main focus of this work will be the so-called unital 
maps for which all bi are zero. Channel $ then maps 
the unit-matrix onto itself and describes a contraction of 
the Bloch sphere without displacement. Such channels 
describe evolutions with decoherence but no relaxation. 
In this unital case, the coefficients obey [20j | . 

| Cl ±c 2 | < |l±c 3 |. (5) 



C. LG correlation functions 

We assume that the time evolution of the system from 
ti to i 2 is described described by the map p — > <&i{p] and 
that from t 2 to £3 by p —> <f> 2 [p] ■ This we assume to hold 
for the measurement of all three correlation functions CV, 
in Eq. ([1]). In particular, the time evolution from t\ to 
t 3 in the correlator C31 is given by the composite p — > 
*i > 2 [^l [p]] ■ This situation is not the most general, as it 
implies the absence of entanglement of the system with 
the bath at time t 2 . This model includes all Markovian 
environments. 

With density matrix p represented by vector \p)), the 
action of each channel map and each measurement su- 
peroperator can be written as matrices acting from the 
left, e.g. ^ ®i\P))- I n tm s representation, the trace 

operation is effected by multiplication from the left with 
the vector = 1(1,0,0,1). F° r generic superoperator 
A, we have Tr{^4[p]} <H> ((l|^4|p)). With measurements 



Sa[p] ^ Qa\p)), this notation allows us to write our ex- 
pressions for the three LGI correlation functions as 

C21 = «l|*aGa*iei|p» 

C 32 = ((l|Q 3 $2Q2$l|p)) 

C31 = ©3*2*1 dip))- (6) 

Since we are not explicitly interested in the unitary 
part of the evolution, we separate each $; as in Eq. ([3|) 
and redefine the measurement operators and initial state 
to remove the unitaries as far as possible. After so do- 
ing, the only unitary part that remains can be expressed 
in terms the single rotation w = u 2 v\. With definition 
of the superoperator, W]/o)) O wpw^, and its inverse, 
W -1 |/°)) w^pw, we obtain 

Cai - ((l|Q 2 5iQi|p)) 

C 32 = ((tlQsW-'S^Q^ilp)) 

C 31 = {(l\Q 3 W- 1 S 2 WS 1 Q 1 \p}}. (7) 

Our task now is to maximise K of Eq. (TTJ) with these 
correlation functions over all measurement angles from 
the parameterisation Eq. This task in too difficult 
to be performed analytically in the most general case. 
We therefore consider some reduced cases of the greatest 
interest. In all cases the maxima are found by setting 
the derivatives of K with respect to the various angles to 
zero and solving. Extensive comparison with numerical 
results was also carried out to ensure that the maxima 
described are global and not local. 

III. SINGLE UNITAL CHANNEL 

We begin with a simple case, in which the second chan- 
nel is purely unitary: S 2 — > 1. In this case, the rotation 
W cancels from all correlators, which then involve only 
the measurements and the Si matrix. With Si unital, 
the three correlation functions in terms of the angles of 
Eq. © read 

Cfei = C3 cos 9 1 cos 9 k 

+ ( ci cos 4>i cos 4>k + c 2 sin (f>i sin (f>k ^ sin 9i sin 9k ; 

C 32 = cos 2 cos 3 + cos (4> 2 - 4> 3 ) s'm9 2 sin0 3 , (8) 

for k = 2,3. These correlation functions, and thus the 
value of K . do not depend on the initial state p (as has 
been observed for the LGI without dephasing |26l]). 

With these correlation functions, K shows local max- 
ima at values K = 1 + [13] ■ The maximal value for 
a given set of c; parameters is therefore 

tfWx = l + §max[c?]. (9) 
Thus, providing that the channel can support some de- 
gree of coherence, i.e. the channel does not map every 
input state onto the maximally-mixed state, then some 
degree of violation of the LGI is possible since Kmix ^ 1 
provided at least one c, is non-zero. The maximization 
over all measurement angles serves to pick out the direc- 
tion in which Cj is largest to obtain the largest violation. 
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A. Example: Mach-Zehnder Interferometer 

The effect of decoherence on the test of LGI with Mach- 
Zehnder geometry [l(], [H| can be formulated within this 
single-channel language. If wc assume that the decoher- 
ence of the two paths occurs exclusively before the second 
measurement, then the decoherence can be modelled by 
setting C3 = 1 and C2 = c±, with the value of C\ depending 
on the strength of the decoherence. 

According to Eq. ©, the maximal violation obtained 
with this set-up should be K max = 3/2 since max [cj] = 1. 
This is the same as in the absence of decoherence. In 
Ref. [ll|, however, decoherence was found to reduce the 
value of K away from this value. The resolution is that 
in Ref. [ll|, the measurement Q2 occurs in a fixed basis, 
namely the basis in which the off-diagonal elements of 
the density matrix are reduced by decoherence. By fixing 
62 = in our work, the three LG correlators read 

C 2 i = cos 6i ; C 32 = cos # 3 ; 

C31 = cos 9i cos #3 + c\ cos (4>i — 4>s) sin d\ sin 6*3 (10) 

which (up to trivial redefinitions) are the same as found 
in Ref. [111 ]. Maximisation of K over the remaining an- 
gles gives A'max = (1 + \ci\ + Ci)/(1 + |ci|), which is 
equivalent that found in Ref. [ll|. This shows that this 
simple single-channel model can reproduce the results of 
a proposed LGI test. 



IV. TWO UNITAL CHANNELS 

We now consider the case where both channel matri- 
ces Si are the same. Physically, this would most likely 
arise when the same decoherence mechanism is active 
throughout the experiment and the times in the LGI are 
equally spaced: ti = t\ + r and £3 = t\ + 2r, as is the 
most-often-studied case. 

As in the foregoing, the specification of unital channels 
makes K independent of the initial state. We will study 
two cases separately here: first the case where the unitary 
component coming from the channels is absent: W = 1 
(which implies implies u = v); and then the case without 
this restriction. 



A. No W rotation 

With W = 1, the maximization of K can be performed 
analytically. Let us denote c = max{c;}, and d the next 
largest of the three ct. Three maxima of K are relevant: 

K^l = \c\{2-\c\), (11) 

which is valid for all parameter values; 

^ = l + C 2 (i--T^), ( 12 ) 
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FIG. 1. The maximised LG quantity if ma x for the two- 
channel unital case with no additional rotations, W = 1, 
plotted against the maximum decoherence parameter c = 
max{ci}. The blue points show AT m ax values calculated for 
a set of randomly generated set Ci parameters. The orange 
dashed line shows the lower bound -ft'max and the red solid line 
shows the upper bound ifmax • Only for a decoherence param- 
eter c > l/y/2 can .Kmax exceed unity such that violations of 
the LGI may occur. 



which is only valid when J^-i < 1; and finally 

<r + c! 

which is only valid for c 2+ c , i i < 1 [28j - Apart from the 
constraints mentioned, the value of A' max is the maximum 
of these three values. 

Figure [TJ shows the value of K mllx obtained for a set 
of randomly generated Ci values, plotted as a function 
of c. In generating the set of Ci values, Eq. ([5]) must be 
obeyed. For |c| < 1/2 the maximal value depends only on 

the value c, with K m&yi = Km&l < 1- Above this value, 
the two further solutions come into play and, for a given 
value of c, K max can take on a spread of values. The 
lower bound of this spread is simply isTmax, the upper 
bound is found by taking ifmlx with c = d which gives 

K^ = \ + c 2 , (14) 

for |c| > 1/2. 

The violation of the LGI has a very clear threshold in 
this case — the upper bound for K mayi only exceeds unity 
if c is greater that Thus, for violations of the LGI to 
occur, it is a necessary condition that 

c = max{ci} > (15) 
v2 

and if dephasing is too strong, violation of the LGI can 
not occur. Eq. (|15p is in itself not sufficient to ensure vio- 
lation, as the lower bound lies within the classical bound 
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FIG. 2. The same as Fig.[T]but with a finite W rotation. In 
obtaining the blue data points, random values of the angles 
6w and <f)w are generated randomly alongside the decoherence 
parameters Ci The upper bound (solid red line), and hence 
the threshold for LGI violations, is the same as in the W = 1 
case. The lower bound is given by the minimum of either 



K 



(2a) (orange dashed) or K^l (pink dashed). INSET: The 
two-unital-channel LG quantity A" m ax with C2 = C3 = and 
4>w = as a function of rotation angle 6w ■ Shown are results 
for values of ci = 0.6, 0.7, 0.8, 0.9, 1, as marked. 



for all values of c. Consideration of the other parameters 
must be made to determine whether there is a violation 
or not. One general remark can be made, however. If 
one or two of the C, are zero, then the maximum value of 

-Kmax is restricted to the solution A'mfi < 1 and there can 
be no violation. This arises from the restrictions placed 
on the Ci values by Eq. ([5]). 



B. With W-rotation 

The addition of the rotation W complicates the max- 
imisation of K significantly, and an expression for the 
complete dependence of K max on all parameters is not 
forthcoming. Nevertheless, some exact statements can 
be made. Figure [2] shows K max for random values of the 
rotation angles 6\y and <pw as well as Ci. For c < 1/2, 
the solution i^max becomes the upper bound for K max , 
and a lower bound of 



K<&> = V2c 



(16) 



develops [29j. For c > 1/2, the upper and lower bounds 
on -Kmax are exactly the same as in the W = 1 case. Im- 
portantly, this leaves the threshold behaviour unaltered: 
a c of greater than 1 / y/2 is required for violation of the 
LGI to be possible. 

An illustrative, exactly-solvable, example of the de- 
pendence of the results on the relative-unitary an- 
gle 0\y obtains when c\ is the only non-zero coeffi- 



by the maximum of |ci| (2| cos6*n/| — |ci| cos2#w) or 
|ci| (2| sin0wr| + |ci| cos2#wr). The behaviour of this 
maximum is illustrated in inset of Fig. [2] and shows a 
strong variation with angle 9w, even ranging from 1 to 
3/2 in the limiting case of c = 1. 



C. Example: Depolasing channel 

Let us consider that the channel in question is a depo- 
larising channel, i. e. unital with c; = c. In this case, we 
hnd 



K — 



|c|(l 



|c|) |c| < 1/2 
• 2 \c\ > 1/2 



(17) 



irrespective of whether W = 1 or not. In this case, K max 
simply follows the upper bound in Fig. [TJ That K miiX is 
independent of W as can be explained by the isotropy of 
S. 



D. Example: Damped Rabi oscillations 

The typical experimental system in which the LGI has 
to-date been investigated is a qubit undergoing damped 
Rabi oscillations. A complete channel description, in- 
cluding unitary part, of this dynamics, reads 



<I> = 



/ 1 
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V - e -2 r *sinm 
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(18) 
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nt ) 



e~ r K to- 



w/2 



cicnt and 



9 lr 



0. In this case, K, 



max is given 



where fi is the Rabi frequency, F the decoherence rate 
and t, the time. This evolution corresponds to a uni- 
tal 5* matrix with c-2 = C3 = e _ 2 rt and c\ 
gether with unitary rotation W with angles (j>w 
and 6w = Qt. 

Maximisation over all measurement angles gives the 
same result as for the depolarising channel, Eq. (|17|) . in- 
dependent of rotation angle Ow- This follows because 
W only acts in a subspace in which the two c-values 
are isotropic. The threshold for violation of the LGI, 
Eq. (fT5)) . therefore holds exactly for this important case 
and, in terms of the original parameters, translates as 
e~ 2 rt > 1/ With the time of maximal violation given 
by nt = tt/6 [J, we obtain a condition on the ratio of 
decoherence rate and Rabi frequency 

£ < -log2« 1.32, (19) 

Si 7T 

for violations of the LGI to be possible. 



V. NON-UNITAL CHANNELS 

The general non-unital case is markedly more complex 
than the foregoing and we will not attempt to give a 
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FIG. 3. -Kmax in the non-unital case. The green crosses 
show values for randomly generated channels with W = 1; 
blue circles, values for random parameters with finite W . The 
solid lines show the bounds without W rotations; the dashed 
lines, the bounds without. For W = 1 violations of the LGI 
occur only when c > 0.544. For arbitrary W, no threshold is 
observed and violations are possible at all values of c 7^ 0. 



comprehensive analysis here. Nevertheless, it is worth 
highlighting a few important aspects that differ from the 
unital case. 



A. Single non-unital channel 

Repeating the single-channel calculation of section IIIII 
with a non-unital S'-matrix, we observe — and this is 
a general feature of the LGI for non-unital evolution — 
that the value of A, and hence A max , does depend on the 
initial state. 

In the single-channel case, however, if the input state 



is the complete mixture p = 1 , it turns out that the LGI 
correlation functions are exactly the same as in Eq. 
This means that for this initial state, the non-unital re- 
sults are independent of the parameters bi . 

In general, though, A depends both on the initial state 
and on the non-unital parameters bi. By way of example, 
consider a "classical channel" [2(| with b\ = 62 = ci — 
C2 = and C3 = 63 = 1/2, with pure state p = 
|V-') = cosA| f) +sinA| J.) as input. Maximisation over 
measurement angles gives A max = j (5 + cos 2 A), which 
oscillates between the classical bound (A max = 1) and 
maximum violation (A max = |) as a function of input- 
state angle A. 

B. Two non-unital channels 

In the two-channel non-unital case, the LG parameter 
A depends on the input state once again. Here we shall 
only consider a completely mixed state as input. 

Figure [3] shows a set of A- values, together with the 
upper and lower bounds, for this non-unital situation. 
Results for both = 11 and with a finite-VF rotation are 
shown. In the W = 1 case, the bounds can be found ana- 
lytically. The upper bound consists of two sections. The 
first consists of the classical bound A' max = 1 , which can 
be obtained by setting e.g. C\ = c 2 = b\ = b 2 = 0; 
C3 = c; 63 = 1 — c and <fi a = 6 a = 0. The sec- 
ond can be found by setting c\ = C2 = c, b\ = 62, 
C3 = &3 = and then determining the value of b\ as fol- 
lows. We consider the action of S on pure state density 
matrix p = n.a with n = (cos A, sin A cos B, sin A sin B). 
The length of the resulting Bloch vector, I = 

y (61 + ccosA) 2 + (bi + csin AcosB) 2 , is maximal for 

input angles (A,B) = (7r/4,0). By setting c + V2b 1 = 1, 
we obtain I = 1 and it is with this condition that max- 
imal value of K is reached for a given value of c. The 
expression for the upper bound so determined is 



2c 



c 2 cos < cos 1 



3c 4 



4(c 2 + c 4 ) 



V-l - 2c 2 + 9c 4 + 28c 6 + 9c 8 + 6c 10 - c 12 
2c(l + c 2 ) 2 



The complete upper bound consists of the first value 
(unity) for |c| < 0.544, and Km&! above this value. The 
lower bound shows an asymmetry with respect to the 
sign of c. For c > 0, the bound is simply Amax- For 
c < 0, it is 

KgV = Max [|c|, -1 + i\c\ - 2c 2 ] . (20) 

With W rotation, the bounds are extended somewhat. 
The lower bound becomes Amax of Eq. (|20[) for all c. The 
upper bound in Fig. [3] could only be found numerically 
This was found by keeping just one Ci non-zero, c\ say, 



and then setting 63 = and 62 = yl — (ci + &1) 2 such 
that the length I = y/ (c\ + bi) 2 + 62 + ^3 is equal to one. 
For a given ci , we have to numerically optimise the value 
of bi that gives the maximum violation. With the above 
set, the maximizing angles are Q\ = n/2, <f> a = cf>w = 0. 
The remaining angles, 82, O3 and 6w, are then found 
numerically 

The most significant feature of Fig. [3] is that whilst 
the violation in the W = 1 case still shows a threshold 
behaviour, with finite W it does not. With W = 1, viola- 
tions are only possible for \c\ > 0.544 whereas violations, 
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even if only small, are possible for all c 7^ in the finite 
IF-rotation case. This enhanced violation arises because 
the relaxation present in the non-unital case can lead to 
an increased purity of the output state in comparison 
with that following from the same channel with hi = 0. 



VI. CONCLUSIONS 

We have used here a maximisation over measurement 
angles to investigate the relationship between the LGI pa- 
rameter K and the non-unitary aspects of the dynamics 
as expressed through the quantum channel parameters Ci 
and hi. The rotation W (relative unitary) between the 
two-channels was also seen to play a role. 

With just a single unital channel, the relation between 
K max and the channel parameters is extremely simple, 
Eq. © , and except for the trivial case in which the chan- 
nel maps every state onto the unit matrix, violations 
of the LGI can always be found. In contrast, the two- 
channel unital case shows a threshold behaviour, such 
that no violations occur if the maximum dephasing pa- 
rameter c is less than l/y/2. This allows us to rule out 
violations of the LGI for a large class of environments. 

The non-unital case is significantly more complex than 
the unital and we have only touched on it here. Im- 
portantly, the LG correlation function and its maximum 
value are no longer independent of the initial state. In 
the most general non-unital case, the threshold behaviour 
disappears, although this re-emerges for restricted classes 
of channel. In the unital case, considering K max as a 
function of c provides a useful way of organising the be- 
haviour of the LGI. For the non-unital case, it is not clear 
whether this is still the case, since the effect on K max of 



the displacements b% are on an equal footing to those 
caused by the contractions Cj. More work is certainly 
needed on this front. 

Our results set upper bounds on the possible values 
of K for a given system. If in experiment the measure- 
ment basis is restricted, or equivalently if the unitary dy- 
namics of the system is not completely controllable, then 
these additional restrictions may mean that the maxi- 
mum value of K obtained in experiment is less than the 
bound -fTmax here. Knowledge of the maximal violation, 
however, may also suggest ways in which the experimen- 
tal measurements could be modified to maximise the vi- 
olation. One further perspective on our results is that 
a measurement of K max allows us to infer certain infor- 
mation about the channcl(s) in question, in particular, 
bounds on c. In this context, the local maxima of K, 
and not just its global maximum, can provide additional 
information, e.g. in the single channel unital case, know- 
ing the local maxima provides all three Ci. Whilst this is 
an interesting perspective on the violations of the LGI, 
it would not represent an efficient way to experimentally 
characterise the channels(s) in question. 

Finally, we mention that the system-environment in- 
teraction considered here is not the most general, as we 
have excluded the possibility of system-bath coherence at 
the point where the second measurement is made. The 
influence of strongly-coupled quantum environments on 
the LGI remains to be seen. 
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